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The usual nuclear recoil energy reconstruction employed by liquid xenon dark matter search 
experiments relies only on the primary scintillation photon signal. Energy reconstruction based 
on both the photon and electron signals yields a more accurate representation of search results. 
For a dark matter particle mass m x ~ 10 GeV, a nuclear recoil from a scattering event is more 
likely to be observed in the lower left corner of the typical search box, rather than near the nuclear 
recoil calibration centroid. In this region of the search box, the actual nuclear recoil energies are 
smaller than the usual energy scale suggests, by about a factor x2. Recent search results from the 
XENON100 experiment are discussed in light of these considerations. 



Introduction.— Liquid xenon detectors are presently at 
the forefront of direct searches for galactic particle dark 
matter. They have placed the most stringent upper lim- 
its on the dark matter — nucleon cross section a n pQ, 
and an additional factor x5 (or more) improvement in 
sensitivity is expected in the coming year [2]. The ex- 
pected signature of a dark matter interaction with target 
nuclei is, in a large number of theoretical models, a low 
energy nuclear recoil [3] . Liquid xenon offers an approxi- 
mate factor x200 discrimination between nuclear recoils 
(as from neutrons, and expected for dark matter), and 
electron recoils (as from electromagnetic background) [4] . 
The basis for the discrimination arises from the fact that 
the partitioning of electronic energy losses into a num- 
ber of photons (n 7 ) and electrons (n e ) is a function of 
the incident particle type, and also of the particle energy 
[5]. Thus search results from liquid xenon detectors are 
generally given in terms of a nuclear recoil energy (E nr ) 
estimator versus the discriminant, y cx n e /n 7 . 

In practice, this has exclusively taken the forms 
E nr cx SI and y cx S2/S1. In these equations, 
SI = fti n 7 and S2 = a 2 n e are the experimentally 
measured number of photoelectrons (from photomulti- 
plier tubes), corresponding to the photon and electron 
signal. Typical values are ol\ ~ ^ and ct2 ~ 20, and 
these quantities are detector-dependent. The purpose 
of this Letter is to show that an energy reconstruction 
based on both n e and n 7 is significantly more accurate 
and useful than the usual E nr cx SI method, especially 
in considering possible detection scenarios. 

Discussion — Nuclear recoil energy reconstruction 
in liquid xenon has historically been defined by 
E nr = (SI) /(Ly£ e ff)(Se/S n ), where C eff is the rela- 
tive scintillation efficiency of nuclear recoils (also referred 
to as the effective Lindhard factor). A monoenergetic 
gamma (usually 122 keV), with a measured light yield 
L y given in photoelectrons/keV, is used as standard 
candle. An electric field applied across the xenon target 
allows the measurement of the electron signal, and it 
also quenches the scintillation signal. The quenching of 
the gamma is S e ~ 0.5, while nuclear recoils are barely 



quenched (S n ~ 0.95). Both L y and S e vary significantly 
with the electric field strength. 
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FIG. 1: (upper panel) Simulated nuclear recoil calibration 
data from a hypothetical detector which is similar to the 
XENON100 detector during its 2009 commissioning run. The 
dashed curves, described in the text, define three sides of the 
signal search box. The inset contours indicate (E nr ). Electron 
recoils from electromagnetic background have (y) = 0. The 
energy scale used in [T] is indicated along the top, and the 
two events reported therein are reproduced for comparison, 
(lower panel) SI acceptance of the hypothetical detector 
(solid), and the XENON100 detector (dashed). 

Dedicated, direct measurements of the relative scin- 
tillation efficiency C e ff are plentiful [6HH]- These mea- 
surements report, for a given E nri the average number of 
recorded photoelectrons (SI). In contrast, dark matter 
search experiments report a measured SI (and S2) corre- 
sponding to each event, and wish to know the most likely 
nuclear recoil energy (E nr ) associated with that event. 
For events near the nuclear recoil band centroid (in ?/), 
these two operations commute. But in many cases, such 
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as the recently reported two events in the XENON100 
search region [1 j, they do not. 

It appears that a nuclear recoil energy reconstruction 
based directly on Lindhard theory is possible in liquid 
xenon [5 j. The method employed in this work is very 
nearly equivalent to that, but framed so as to make an 
explicit connection with the historical method, and with 
the experimentally measured quantities. 

Figure [I] shows the results of a simulation of the nu- 
clear recoil response of a hypothetical liquid xenon detec- 
tor. The results are plotted in the customary variables, 
though other choices are possible [9], and may be useful 
considering the relative size of oq and ol^. The simulation 
method is described in [10], and reproduces the relevant 
binomial and Poisson statistical processes. It has been 
shown to provide a hi-fidelity reproduction of XENON10 
nuclear recoil calibration data. The detector-specific de- 
tails were obtained from [TTJ [12] . The hypothetical de- 
tector is therefore expected to exhibit a response simi- 
lar to the XENON 100 detector during its initial phase 
of operation. Uncertainties are discussed in a separate 
section. In this work, we took as inputs the nuclear re- 
coil centroid reported in Fig. 3 of [11], and the central 
C e ff curve from [13] (shown in Fig. [2j lower panel, solid 
curve). The simulated data are shown after subtracting 
the electron recoil centroid (/aer), which was also taken 
from [TT] . 



(solid curve). This does not guarantee that either £ e // 
or Q y are correct in absolute terms, but rather, as drawn, 
are self-consistent with nuclear recoil band data. We note 
that the lower (dashed) Q y curve is very consistent with 
the NEST model |E]. 

In addition to the simulated nuclear recoil data, Fig. [I] 
shows the centroid /anr (dashed, green) and Hnr — 3&nr 
(dashed, black). Also shown are the "software" SI > 3 
(dashed) and S2 > 150 (stippled) thresholds, as in pQ. 
The dashed curves define three walls of the dark matter 
search box for our hypothetical detector. In a similar 
(actual) search box, the XENON100 Collaboration re- 
cently reported the observation of two events, which are 
reproduced here. The stated event energies are 7.1 (cir- 
cle) and 7.8 keV (star), and this appears very reasonable 
according to the E nr ex SI scale given along the top axis. 
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FIG. 2: Nuclear recoil energy calibration data for liquid 
xenon. The simulation assumes the solid curves. Events with 
E nr < 2 keV were not simulated. C e ff data reproduced from 
[6] (stars), [7] (diamonds) and [8] (squares). Additional Q y 
data from [15] . 

The simulation has only a single free parameter, the 
ionization yield Q y = S2/E nr . It was allowed to float 
until the simulation band centroid matched the nuclear 
recoil centroid from data [11 . The agreement is very 
good, within la of the statistical uncertainty on the 
mean, above SI = 3. Below SI = 3, the agreement is 
within 2<j. The Q y curve so obtained is shown in Fig. [2] 



FIG. 3: The probability to observe a fluctuation equal to or 
greater than the two events shown in Fig. [I] as a function of 
nuclear recoil energy. The abcissa is the actual simulated en- 
ergy (not the contours in Fig. [T]), and the markers correspond 
to the two events. 

We used the simulated events to find contours of (E nr ), 
as shown in Fig. [I] (solid curves, with corresponding en- 
ergy in keV). Because a downward fluctuation in SI is 
accompanied by an upward fluctuation in y, the contours 
follow the S2 expectation value for each energy. This is 
significantly different from the cartesian expectation im- 
plied by E nr oc SI (notice that this scale is most correct 
near the calibration centroid, near y « -0.4). Figure [3] 
shows the probability that a nuclear recoil of energy E nr 
resulted in either of the two observed events. Only sim- 
ulated events which produced a measurable SI and S2 
signal were considered. 

The discussion up to this point has assumed a nuclear 
recoil spectrum corresponding to an americium-berylium 
neutron source, as is frequently used to calibrate the nu- 
clear recoil response of liquid xenon detectors. In Fig. 
[4] we show the expected distributions of events for sev- 
eral dark matter masses m x . The spectral distributions 
were calculated assuming the same astrophysical param- 
eters described in pQ. This clearly shows, particularly for 
m x < 10 GeV, that the y coordinate also carries spectral 
information. The point is perhaps obvious from the def- 
inition of but it has been neglected in previous work. 
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FIG. 4: The expected distribution of nuclear recoil events, for several values of m x . In order to clearly show the distribution, 



Uncertainties.— The ordinary statistical processes mod- 
eled by the simulation lead to non-Gaussian tails in the 
distribution of y, particularly for SI < 10. The extent 
of the tails are roughly indicated by the (E nr ) contours 
in Fig. [I] In Fig. [5] (upper panel), we indicate the ap- 
proximate range of Q y (dashed curves) which adequately 
reproduce the nuclear recoil band data centroid in [TT] . 
given C e ff as in the lower panel (solid curve). The 
lower dashed curve is the model prediction employed in 
[21]; Q y below this do not appear reasonable. It can be 
seen that the uncertainty in Q y is most significant for 
E nr < 10 keV. This is due primarily to the non- Gaussian 
tails in y, and how the mean value of the distribution is 
determined (our method may be slightly different from 
what was used in [H]). As a result, one may expect a 
systematic shift in the (E nr ) contours in Fig. [I] of as 
much as Ay =to!o6 for E nr < 10 keV. 

The SI acceptance of our hypothetical detector re- 
sponse is different from pQ, as shown in Fig. [I] (lower 
panel). The statistical methods we employ make it easy 
to adjust the predicted SI acceptance to match that of 
PQ. This has a systematic effect on the best-fit Q y , as 
pointed out in [16], and leads to a maximum displace- 
ment Ay = +0.02 of the (E nr ) contours. This effect 
causes the largest shift in the range 3 — 5 keV. 

Of course, there is also uncertainty in C e ff itself. In 
this work, we have taken as a prior the C e ff favored 
by the XENON100 Collaboration. An equally plausible 
choice would be [6] (shown in Fig. [2j stars). A 20% 
smaller C e f / would require a systematic decrease in Q y of 
about 8 — 12%. This would lead to an essentially uniform 
—5% systematic shift in the (E nr ) contours, as can be 
verified analytically. 

It is also notable that the SI response reported in [1 
has improved by ~ 4% and the S2 response has improved 
by ~ 14%, relative to [11 (and hence relative to our hy- 
pothetical detector). This leads to a uniform +4% sys- 
tematic shift in the (E nr ) contours. We mention in pass- 
ing that the width of the simulated band (in y) appears in 
good agreement with previously reported results [TTJ [13] , 



but slightly narrower than [I. In spite of these uncer- 
tainties, the fundamental profile of the (E nr ) contours 
remains as shown in Fig. [I] 

Conclusions.— Full consideration of the energy informa- 
tion carried by both the SI and S2 signals indicates 
that the energies of the two events observed by the 
XENON100 detector would be reconstructed at 2.9 ± 
0.5 keV and 3.6 zb 0.6 keV in our hypothetical detector, 
subject to a systematic uncertainty of l^o keV. This is 
about a factor x2 smaller than what one obtains from 
E nr oc SI. Interestingly, the approximate location 
of these events (near the lower left corner of a typical 
search box, rather than near the nuclear recoil centroid) 
is what one would expect for elastic scattering of low- 
mass (m x < 10 GeV) dark matter. This observation 
could be important from a phenomenological perspective 

(e.g. [mm)- 

It is evident from Fig. [4] is that the acceptance of the 
search box (as a function of SI) depends on m x . For 
the four example m x values, the fraction of events below 
the calibration centroid in the range 3 < SI < 10 are 
1.00, 0.99, 0.56 and 0.42. This is a result of the shape of 
the calculated dark matter energy spectrum. The frac- 
tion is about 0.53 for the calibration data, slightly above 
0.50 due to the non-Gaussian tails. Energy reconstruc- 
tion based on E nr oc SI leads to the assumption that 
the acceptance of the search box is always given by the 
calibration data. Relative to previously reported results, 
this should tend to strengthen the sensitivity to parti- 
cle masses m x < 35 GeV, and weaken the sensitivity to 
larger masses. From an experimental point of view, it is 
interesting that for m x < 10 GeV, the electromagnetic 
background population with (y) = is essentially irrele- 
vant. 

In Fig. [5] we show values of (m x ,a n ), in which out- 
comes similar to that recently observed by XENON 100 
appear at least 10% likely (crosses) in our hypothetical 
detector. The SI acceptance is shown in Fig. [I] (solid 
curve), and a 34 kg x 224.6 day exposure was assumed. 
Our definition of "similar" is two events observed in the 
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FIG. 5: The region of parameter space in which an outcome 
similar to the two events shown in Fig. [I] is at least 10% likely 
(crosses), and the 90% CL exclusion limits from XENON100 
[1] (thick curve) and XENON10 [21] (thin curve). (Inset) the 
result of a simulated experiment, for a n = 3 x 10 -43 cm 2 and 
m x = 7 GeV. The two simulated events are shown as large 
circles; regions (a) and (b) are described in the text; markings 
and axes are as described in Fig. [I] 

search box below the (E nr ) = 5 keV contour, and no 
events above it. The search box is shown in Fig. [I] 
and Fig. [5] (inset), bounded by dashed lines. Addition- 
ally, to account for the uncertainties discussed above, we 
allow (a) the possibility that one of the two events ap- 
pears in the region below the search box, but above the 
S2 threshold, and (b) the possibility that one of the two 
events falls in the region Sl< 3. These regions are indi- 
cated in Fig. [5] (inset). Uncertainties would propagate 
into Fig. [5] as follows: a larger SI acceptance (e.g. Fig. 
[l] dashed curve) would tend to push the (m x ,<j n ) region 
to smaller <j n . A larger Q y would mean a smaller energy 
for events near the lower left corner of the search box, 
and would therefore tend to push the (m x ,cr n ) region to 
smaller m x . 

We do not suggest that these two events observed in 
P] are due to the elastic scattering of dark matter; the 
background hypothesis is of a similar likelihood, and thus 
more compelling. However, we have shown that in con- 
sidering detection scenarios, significant additional infor- 
mation is gained from an energy reconstruction based 
on both n e and n 7 . Specifically, while a result simi- 
lar to [1 is compatible with low-mass dark matter, it 
is highly unlikely to have arisen from dark matter with 
m x > 10 GeV. 
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